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Abstract 



Often a number of precipitation processes in steels and alloys occur si- 
multaneously (they have the same origin in time) albeit at different rates. 
Consequently isothermic transformations are accompanied, in this case, by 
the occurrence of second processes of precipitation that influence the changes 
of the macroscopic parameter chosen for the kinetic study while first processes 
are occurring. In this context a set-up is presented that addresses the issue of 
characterization of nucleation and growth mechanisms of the individual pro- 
cesses. Soft and hard interferences of processes are considered. The present 
approach is based on the Avrami model and allows computing the kinetic 
parameters (n, k) of the different precipitation processes in isothermic phase 
reactions. 
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I. INTRODUCTION 



The Johnson-Mehl-Avrami(JMA) equation has been widely used to study the kinetics 
of the nucleation and growth of a second phase from an initial parent phase [1-5]. This 
kinetic relation needs assumptions about the nucleation and growth rates of the new phase 
nuclei [6-10]. By evaluating the JMA equation, an index n is obtained. Knowledge of n, 
together with the empirical value of the activation energy and some other features of the 
transformation, yields a much better understanding of the reaction mechanism. 

In practice, a macroscopic parameter P(t) is commonly used for the kinetic study of 
the isothermic reactions [10]. By simplicity, it is assumed that P(t) depends linearly with 
the solute concentration c(t) that participates in the reaction. Variations of P(t) during the 
transformation are only associated with a type of change; either by entry or exit of the solute 
in the parent phase or by formation /dissolution of a second phase. The phase fraction that 
is transformed in time can then be written as: 



Often, in practical situations as the low-alloy steels martensite tempering and many other 
isothermal phase transformations in steels, precipitation processes do not occur isolately 
[11-13]. Instead, first and second processes develop simultaneously (the same origin in 
time) albeit at different rates. Then, the macroscopic experimental parameter P(t) contains 
information of both processes (overlapping processes). In this sense, some authors [14] 
have remarked evident difficulties while trying to set apart the effects of each one of these 
processes by knowledge of only the macroscopic parameter P(t). 

Most of these works described above have involved fitting of an Avrami type equation 
to the transformations curves for a specific alloy at specific temperatures. These have con- 
sidered only one precipitation reaction neglecting the influence of the other precipitation 
processes from the experimental data. Therefore, the variation of the fraction transformed 
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over time is frequently found to deviate from a simple sigmoidal relationship ( like JMA be- 
havior), because the curve vs. hit actually represents more than one reaction occurring 
at the same time. 

There have been a number of attempts to deconvolute such master curves into compo- 
nents due to individual processes [15-18], but it is difficult to see whether some of these 
fitting parameters have physical meaning. 

The aim of the present paper is, precisely, to outline a set-up based on the Avrami model, 
that allows obtaining the kinetic parameters (n, k) for the simultaneous precipitation pro- 
cesses with interference by the only knowledge of a macroscopic parameter P(t) chosen for 
the kinetic study. This new set-up will be helpful for the interpretation of the experimental 
data. It does not fall into the category of the predictive models of references [19-21]. 



II. KINETIC THEORY OF TWO OVERLAPPING PROCESSES. 

Assuming that JMA equation holds for individual processes is equivalent to considering 
only hard impingement between particles that participate in each process. In other words, 
the soft impingement for a given type of particles in the same spatial region is not considered. 

Let us take two processes that rival for the same solute at different rates. Suppose both 
obey JMA equations. Then, the real fraction transformed that is determined by the changes 
of the experimental parameter obeys a relation that is similar to Eq. (1): 



where t en d\ is the time at the end of the first process. In Eq. (2) we chose t en d\ instead of 
tend2 because we take into account the effect of the first process during the interval t < t en d\ 
when there is appreciable influence of the second process on the experimental parameter. 
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When the existence of a second process is known, it is used to fix the end of the reaction at 
the time when the process that causes major influence in the parameter P(t) ends. In the 
case when this last is the second process, it is used to neglect variations in P(t) due to the 
first one. Obviously it is incorrect. 

It is interesting finding the relation between £ r (t) and the fractions £i(t) and £ 2 (£) of the 
elementary processes by knowledge of AP(t) only. 

Let us consider that variations of P(t) associated with each process are independent, i.e., 



Therefore, we neglect the reciprocal interference of both processes, i.e., the overlapping 
of the drainage zones of the chemically different precipitates that correspond to processes 1 
and 2. The bigger the differences of the rates at which the individual processes take place, 
the higher the reliability of this approximation. Neglecting the correlation effects in Eq.(3) 
is less significant than neglecting soft impingement for precipitates of a same kind since, 
phase transformations in both processes, often, take place at different geometrical places. 
In case that the rate of one of the processes is close to that of the other process, rivalry for 
the same solute leads to the situation when a phase alone rivals for the same solute, in other 
words, a situation with soft impingement among particles of a given kind. But we recall that 
the present set-up, soft impingement are neglected, i.e., JMA-type expressions are satisfied. 

Therefore, for each elementary process we have: 



AP(t) = AP 1 (t) + AP(t 2 ). 



(3) 



AP 1 {t) 



(4) 



APi(t en di) 



and 



AP 2 (t) 



(5) 



AP 2 (f erK 2 2 ) ' 
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respectively so, 



AP 1 (t) + AP 2 (t) 



(6) 



APi(t en dl) + AP 2 (t end i) 



If we take into account Eqs.(4) and (5), then Eq.(6) can be written as, 




(7) 



where 



AP2(t en d2) 
A Pi (tendl) 



(8) 



The magnitude ^(tendi) measures the degree of overlapping (soft or hard) of both pro- 
cesses. Depending on the choice of P(t), the parameter a p may be either negative or positive. 
When P(t) moves in contrary senses for both processes, a p < but when P(t) moves in the 
same sense, a p > 0. 

We remark that £ r (t) depends on the parameter chosen for the kinetic treatment of the 
processes. Recall that a p takes different values for the different macroscopic parameters P(t) 
one chooses. 

A. Computation of £ r (i) with the help of the isothermic dilatometric records. 

Dilatometric records are often used for the characterization of the transformation kinetics 
of a significant group of steels [3,22-24]. The experimental determination of £ r (i) would be 
possible if we would be able to measure the parameter P(t)=l(t) from the beginning of the 
phase transformations (in an isothermic regime). However, the sample needs some time 



5 



for reaching the temperature of the isothermic treatment (for example, the temperature of 
isothermic tempering of martensite structures). During this small time interval, the sample 
suffers a thermal treatment, so the beginning of the transformations may be prior to that of 
the temperature stabilization corresponding to the isothermic treatment and therefore, the 
beginning of the precipitation precesses is unknown. 

Without loss of generality, the figure 1 shows a scheme corresponding to a dilatometric 
record belonging to a low alloy steel martensite tempering. In the interval < t < t (fig. 
1), a dilatation of the sample occurs that is originated by the change in the temperature 
from a room temperature up to the temperature of the isothermic treatment. Simultane- 
ously, a contraction of the sample takes place in this interval that is a consequence of a first 
process in the chain of transformations during the heating of the sample. Consequently, it 
is necessary to begin the study not from the origin of the record given by the measuring 
equipment but from the moment of time when the isothermic regime is reached. It is given 
by l' (0) in fig.l. If we take /'(O) to be the initial length of the sample (/'(0) = Z(£ ))then we 
may compute the fraction transformed through: 



We remark that £(t') is not the real fraction transformed £>(£). This last magnitude 
may be written, following Eq.(2) with P(t)=l(t) as: 



Since £ r (t)(Eq.(10)) can not be obtained directly from experiment (we do not known 
the origin (in time) of the transformations), it is necessary to find a relation between the 
fraction transformed £(t'), that is experimentally measurable, and £ r (i) . It can be shown 
that: 
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where: 



z\t) = z\t-t ) = at)- (12) 

From Eq.(ll) we obtain: 



d£(t') _ d£(t') 



dint' (1 - £ r (t ))d\nt' 



since: 



Ut) = Ut +t ) = C(t). (14) 

If we take the derivative of Eq.(7) and considering that £ r (t) = £ r (t') and £(t') = £(t) 
the following expression is obtained: 



but 



dint' 1 + a^ 2 (t' endi y dint' 1 dint' U { } 



dint' ~ t' + t dint ■ [ib) 

Let £ (t) be the fraction transformed of an elementary process obeying a JMA expression, 
£(f) = 1 - exp - (kt) n , then: 



-^^dnVt =nOI JnVt =nZ ^ (17) 
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where 



d(lnln-r-r) , 

z(0 = [ ^ T = -(! - «i - O- (is) 

Therefore, such as £i(t) and £ 2 (£) both obey JMA-type expressions, we obtain: 



yielding that 



^=^(6) ,(i=l,2), (19) 



^rn^ = -r^^Zfe), (20) 
dint' t'+to 



where the condition £'(t') = has been taken into account. Substitution of Eqs.(20) 
and (15) into Eq.(13) yields a general expression correlating £(t') with the fractions of the 
individual processes that obey JMA-type equations: 



(i - Uto))(i + «^(W))^y = j^t M^ + «^(6)]- (21) 

The behavior of the function Z(£) vs. £ is shown in fig. 2. This function is not sym- 
metric respecting its maximum at £ = 0.632. It grows when £ increases and when £ — ► 1 it 
decreases very rapidly. However Z(£) is nearly constant for £ lying in the neighborhood of 
the point where this function is a maximum. In other words, we have an interval where £ 
takes values (far from the boundary points £ = and £ = 1) for which Z(£) depends weakly 
on £. This is a focal point of our set-up. If we allow time intervals where and &{t) lie 
far enough from the boundary points £ = and £ = 1, then we may consider that Z(£i) and 
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Z(^ 2 ) are nearly constants. The former requirement implies, also, that t << t' so it may 
be neglected when compared with t' . Consequently Eq.(21) may be simplified: 

(1 - er(*o))(l + "i&C*))^ = m^i) + wZfa). (22) 
Two cases are of particular interest and we shall analyze them in what follows. 

1. Soft overlapping of the processes: 

A second process occurs, but it manifests itself only weakly during the interval t < t < 
tendi (0 < t' < t' endl ), i.e., £2 (^endi) is small so, the second process perturbs the first one. In 
this case, since ^(tendi) is small, we may assume that Z(^ 2 ) ~ for to < t < £ e ndi- 
Therefore, Eq. (22) can be written as: 

(1 - *r(<b))(l + «^(W))f^ = nxZ^i), (23) 



or 



^) = nZ(^), with n = - , , w ni -m --. (24) 

dhiif yqih (l-^(to))(l + «^(C d i)) 

In this case when we have a weak perturbation of the macroscopic parameter by a second 
process in the interval t < t a < t < t b < t endl (or < t' a < t' < t' b < t' endl ), where t a and t b 
(or t' a and t' b )&re far from the boundary points, (or £(t')) is always close to £i(t) and, 
consequently, Z(£i) ~ Z(£) without taking into account the stability region. Therefore, the 
following equation takes place (see Eq.(24)): 
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^r = nZ(0. (25) 

In this time interval there exists a linear dependence of Inlnj^r vs Int' so, the slope of 
this dependence is, precisely, n. In other words, in the time interval where the best fit of 
the Inlrij^p- vs Int' is manifest, the experimental data satisfies a nearly JMA behavior. 

Analogously, for times t en d\ < t < t enc i2, (t' en di < t < t' end2 ), when the first process end 
U P (£1 — 1 an d Z(£i) = 0) we have: 



(1 - £r(<b))(l + = «^ 2 Z(6), (26) 



or 



-77— r = nZ{i 2 ), with n = - , -7—, rr- 27 

d^t (l-CrWJll + aiWendl)) 

Since Z(£i) = 0, the second process develops alone in this time interval, therefore 
Z(£ 2 ) — Z{£), where £'(£') decreases while £ 2 (t) increases for the studied case (cti < 0). 
Then, for a certain time interval in the neighborhood of the maximum of Z(£) the points 
(Z(£'), £') and (Z(^ 2 ), £2) are symmetrically located respect to this maximum. This as- 
sumption enables one finding a time interval where Inln-^r vs Int' is nearly linear. This 
means that the experimental fraction transformed follows a nearly JMA behavior. Eq.(27) 
can then be written as, 



f[£ (28) 



and picking a linear arrangement allows computing n by finding the slope of this curve. 
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Relationships (23), (24), (26) and (27) yield a relationship between the a\ parameter, 
the JMA indexes for each individual process and the values n and n found from linear ad- 
justment of the experimental data: 



ai = — . 29 

n 2 n 



2. Hard overlapping of the processes 

In this case a second process is strongly manifested during the time when the fist one is 
occurring. We start from Eq.(22) for the time interval t < t < t endl . As already remarked, 
it is necessary to pick a time interval such that £ 2 > 0.25, i.e., Z(£ 2 ) is nearly constant 
and, consequently, £i(£) and £'(t') are below 0.9. This means taking t < t a < t from the 
left and t < tp < t en d\ from the right. Therefore for times in the interval [t Q ,^][or t' a ,t'p], 
Z(£i) ~ Z(£ 2 ) ~ Z{£) ~ const and the Eq.(22) simplifies to: 



d£tf) 7(c \ + ain 2 

— — t = nZ(t , with n = — — — --— -j—, tt, 30 

dint' (l-er(*o))(l + a,e 2 (W)) 



Yielding that, in this time interval, the experimental data gets linearized, i.e., Inlrij^p-vs Int' 
is nearly linear. The case when the first process ends up {t endl < t < t end2 ), is just the case 
of soft overlapping already studied (Eqs.(26) and (27)). The relationship between ai, the 
indexes for individual processes {n\ and n 2 ) and those obtained from the linear adjustment 
of the experimental data (n, n) for this case of overlapping can be established with the help 
of Eqs. (27) and (30): 



a t = — — . (31) 

n 2 [n — n) 
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B. Testing the kinetics of the individual processes 



Testing the kinetics of a second process allows one obtaining £ 2 (*) o r ^(0 an d therefore 
the parameters n-i and hi can be computed with the help of the Avrami model. 
Therefore, for t > t endl (t' > t' endl ) , £i(t) = 1, yielding that: 



1 + ai&itendl) 



[1 + a&(tendlMr(t) = 1 + (33) 

Besides, from Eq.(ll) it can be obtained: 

er(0=e'(*')[l-er(<b)]+^r(<b). (34) 

Substitution of £ r (i) from Eq.(34) in Eq.(33) allows one obtaining: 



6W=£ 2 (0 = f K'(0-1] + 6(W), (35) 



where Eq.(27) has been taken into account. The fact is that ^(tendi) is unknown. However 
if we set t = t en d2 in Eq.(35) we obtain: 



6(W) = f [l-^L d2 )] + l, (36) 



so, finally, Eq.(35) can be written as: 



6(*) = &*') = ^ K'(0-e'(C d2 )] + i, (37) 

12 



this expression allows finding £ 2 (f) (or £ 2 (0) by the knowledge of £'(£') from experiment. It 
is just our working expression. 

For computation of £i(t) (t < t < t endl ), from Eq.(7) we have: 

6(f) = £(*') = &•(*)[! + a,e 2 (W)] " (38) 
If we put Eq.(34) into Eq.(38) we obtain the working formula for this case: 

= £(*') = £'(*')[! - Cr(*o)][l + a,e 2 (W)] + Cr(*o)[l + ai6(*endl)] - ( 39 ) 

In what follows we shall study the two possible overlapping situations of the processes 
by the analysis of Eq.(39). 

a) - Soft overlapping. 

For times in the interval (t < t a < t < t b < t endl ) (or < t' a < t' < t' b < t' endl ), and with 
the help of Eq.(24), the relationship (39) may be simplified: 

ei(*) = 6(*') = i + ^(*')-i], (40) 

where we assumed that £ 2 (fendi) ~ 62(f) ~ 0. 

b) -Hard overlapping. 

For the time interval (t < t a < t < tp < t end i), with the help of expression (30), the 
relationship (39) can be written now as: 

£i(t) = £(*') = i + ^^'(0 - 1] + ^7^y[6(W) - Ut)l (41) 
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We recall that, for a time interval (t' endl < t' < t' < t' s < t' end2 ), the second process de- 
velops alone. Therefore, lineal independence of the exponential functions leads to n and n 2 
being well correlated. This can be shown by taking all the functional forms satisfying the 
exponential behaviour £'(t') vs. hit' in the interval [iy,^], in the case of interest; a\ < 0, 
(fig. (1)), and satisfying, at the same time Eqs.(28) and (37) in the forms: 

( J ) fj$=*Z(0, with n<0, 



and 



(in m = i&(0-i] + e'(C 2 )- 

The functions describing the experimental behavior £'(t') vs. hit' for t' £ [t'~,t' s ] are: 



f'(f') = a - exp(-kt') n with n < 0, (42) 



and 



f'(i') = b + exp(-kt') n with n > 0. (43) 
Both satisfy expression (I). For the functional form (42) we have, (see Eq.(18)): 



=n(H') n exp(-fct') n = nZ{£), with n < 0. (44) 

For (43): 



^Ml = -n(kt ) n exp(-H') n = -nZ{£) = nZ(£), with n = -n. (45) 
d In. i 
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It is worthwhile noting that only the functional form (43) satisfies condition II with 
n = n 2 . Therefore, since n = —n this yields that n 2 = —n . 

Once the relation between n and n 2 is found, Eq.(41) simplifies: 

= £(*') = 1 + " - eM-htn (46) 

lb lb 

This equation allows obtaining the kinetic parameters k\ and rt\ of the first process in 
the case of hard overlapping of processes for i G [t' a , t'p\. 

With the further purpose of applying the kinetic theory of overlapping processes, we 
compute n as the slope of the linear arrangement In In 1 _^ t /^ v s In t' for the best fitting 
interval t' G [t' 7 ,t' s ] where Z(g) w Z(£ 2 )- 

Functions £ 2 (t') are generated from Eq. (37) for the same interval by means of an iterative 
procedure with an arbitrary initial value n 2 . A data £ 2 (t') is generated for t' G [t' y ,t' 5 ] and 
by further linear adjustment a new n 2 value is computed. The procedure is repeated until 
a desired accuracy is reached. This way the empiric parameters n 2 and k 2 are generated. 
Hence, £' 2 (t' endl ) can be evaluated following a JMA expression and the kind of overlapping 
of processes can be determined. 

Once the kind of overlapping is known the correct expression (40) or (46) is chosen, in 
order to obtain the kinetic parameters of the first process. 

A similar procedure enables one obtaining n but taking the time interval t' G [t' a ,t' b ] or 
t' G according to the overlapping type. 

Once n, n, n 2 , k 2 and £,' 2 (t' endl ) are known, a first data ^[(t') may be computed according 
to whether kind of overlapping Eqs. (40) or (46) one have for an arbitrary initial value n\ in 
the above fitting interval t' . The data generated is linearly adjusted again to obtain a new 
value n x . The procedure is repeated until the desired accuracy is reached while obtaining 
n\ and k\. 
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III. SUMMARY 



The set-up outlined in the present paper opens up new possibilities for the study of 
isothermic phase transformations in a variety of alloy systems. We require only the knowl- 
edge of just some macroscopic parameter that is sensible to changes in the fraction trans- 
formed. The kinetic parameters n and k can then be computed following a very simple 
procedure. This allows, in turn, obtaining a more complete information about the elemen- 
tary mechanisms of nucleation and growth of chemically different nuclei (processes 1 and 2) 
that are competing for the same solute. 

A key point of this set-up is that computing procedure works only in the interval where 
the function Z(£) is nearly constant, leading to the experimental fraction transformed having 
behavior nearly to JMA type expression. 

In a forthcoming paper we apply the set-up presented here in the particular case of the 
martensite tempering of a Mn low-alloy steel. 
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Figure Captions 

Figure 1. Scheme of a dilatometric curve corresponding to the isothermic tempering of 
a low-alloy steel. 



Figure 2. Functional dependence of Z(£) vs. £. 
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